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Abstract. Two different aspects of the description of three- and four- nucleoli 
systems are addressed. The use of bound state like wave functions to describe 
scattering states in N — d collisions at low energies and the effects of some of 
the widely used three-nucleon force models in selected polarization observables 
in the three- and four-nucleon systems are discussed. 

1 Introduction 

Detailed studies in the three- and four-nucleon systems gives valuable infor- 
mation of the underlying nuclear interaction. These two systems have three 
bound states, 3 H, 3 He and 4 He, therefore much of the efforts have been done 
in the study of continuum states. Although a reasonable agreement with the 
available experimental data is obtained in the description of the differential 
cross section in the low energy region, discrepancies can be observed in some 
polarization observables [TJ [2] . Related to this, the analysis of the effects of the 
three-nucleon forces are of crucial importance. Recently a critical comparison 
of different models widely used in the literature has been performed [3]. 

A different aspect of the problem regards the methods used to describe 
continuum states in few-nucleon systems. In the A = 3, 4 systems well es- 
tablished methods to treat both, bound and scattering states, are the so- 
lution of the Faddeev equations (A = 3) or Faddeev-Yakubovsky equations 
(^4 = 4) in configuration or momentum space and the Hyperspherical Har- 
monic (HH) expansion in conjunction with the Kohn Variational Principle 
(KVP). These methods have proven to be of great accuracy and they have 
been tested through different benchmarks [4j [5]. On the other hand, other 
methods are presently used to describe bound states: for example the Green 
Function Montecarlo (GFMC) and No Core Shell Model (NCSM) methods 
have been used in nuclei up to A = 10 and A = 12 respectively The pos- 

sibility of employing bound state techniques to describe scattering states has 
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always attracted particular attention. Recently continuum-discretized states 
obtained from the stochastic variational method have been used to study a+n 
scattering [8]. In a different approach continuum states have been obtained 
using bound state like wave functions [9]. 

In the present paper we discuss the description oi N — d scattering states 
using bound state like wave functions and we briefly show three-body force 
effects in selected polarization observables in A = 3, 4. 

2 Continuum states from bound state like wave functions 

Following Refs. [91 [10], it was shown that a second order estimate of the 
scattering matrix R at a collision energy E (below the breakup threshold) 
results 

B fj = ~ < w i\ H ~ E \Fj > I R 2 nd = A~ l B 2nd (1) 
Aij = <V l \H-E\g j > J ' w 

The eigenvalues of R 2 " d are second order estimates of the phase shifts and 
the indeces (i, j) indicate the different asymptotic configurations accessible at 
the specific energy under consideration. In particular, for the three-nucleon 
system, J-j and Qj are the channel wave functions describing the possible rel- 
ative states of the deuteron and the incident nucleon. For a given J n state the 
different channels are labelled by the relative angular momentum L between 
the deuteron and the incoming nucleon coupled to the total spin S = 1/2 or 
3/2 obtained coupling the spin = 1 of the deuteron the the sin s = 1/2 of 
the incoming nucleon. Specifically j = L,S,J and the channel functions are 



Flsj = J2 



E 

l a =0,2 



w la (xi)F L ( yi ) \ [Y la (x. 



Qlsj = 



v L Jo * J J z 



l a =0,2 



The sum on i runs over the three cyclic permutations of the jacobi coordinates, 
Xi,yi are their moduli, Xi,y~i their directions and u>i(xi) the l-wave deuteron 
wave function. The functions Fi(yi) and G^Vi) are the regular and irregular 
solutions of the N — d Schrodinger equation outside the interaction region. 
The irregular solution has been opportunely regularized at the origin as 



G L (y) = (l-e-^) L+1 G L (y) 



(4) 



where r^d = (\/3/2) y is the nucleon-deuteron separation and 7 a parameter 
that is fixed requiring that Gl(v) = Gi(y) asymptotically. Moreover, Fl,Gl 
are the regular and irregular Bessel functions or Coulomb functions in the 
case of n — d or p — d scattering, respectively. 

The relations given in Eq.([T]) are derived from the KVP, formulating it in 
terms of integral relations depending on the internal structure of the wave 
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function ^ . In fact (H — E)J-j and (H — E)Qj go to zero in the asymptotic 
region since J~j,Gj are the solutions of (H — E)J r j,Gj = in that limit. 
Therefore, in Eq.(pQ), it would be possible to use trial wave functions W'i that 
are solutions of (H — E)^ = in the interaction region but do not have 
the physical asymptotic behavior indicated in Eq.(|3]). In particular, it would 
be possible to use the bound state like wave functions which are solutions of 
(H — E n )'I r ( n * > = in the interaction region at particular values of the energy 
E n . To explore this possibility, let us define a complete square integrable basis 
| J 71 ", a > to expand a bound state like wave function corresponding to a state 
having total angular momentum and parity J 71 ", 

tf» =Y J A n a \J K ,a> (5) 

a 

The index a indicates all the quantum numbers necessary to define the state 
and the linear coefficients of the expansion can be obtained from the following 
generalized eigenvalue problem 

K' < JW > a \ H - E n\J*, a' >= . (6) 

a' 

For example, considering the state J 71 " = l/2 + of the three-nucleon system, 
the lowest eigenvalue after the diagonalization procedure corresponds to the 
three-nucleon bound state energy of 3 H (T z = —1/2) or 3 He (T z = 1/2). How- 
ever, as shown in Ref. [9], more negatives eigenvalues could appear verifying 
\E n \ < \Ed\, with Ed the deuteron binding energy. The corresponding eigen- 
vectors <Z^ n ) approximately describe a scattering process at the center of mass 
energy E® = E n — E^, though asymptotically they go to zero. Considering 
other J w states, the diagonalization procedure will not produce bound states 
since, in the three-nucleon system, a bound state exists only in the J n = l/2 + 
state. However negative eigenvalues could appear, verifying \E n \ < \Ea\- As 
in the previous case, the corresponding eigenvectors approximately describe 
N—d scattering states, though asymptotically they go to zero. The eigenvalues 
E n are embedded in the continuum spectrum of H which starts at En- Ac- 
cordingly, increasing the dimension of the basis the number of them increases. 
We can consider these states approximate solutions of (H — E n )\p( n ' = in the 
interaction region and use them as inputs in the integral relation to compute 
second order estimate of the phase-shifts. As an example, results for scat- 
tering J = l/2 + ,3/2 + states are given in Fig. [1] using the s-wave MT I-III 
nucleon-nucleon interaction [llj. The n — d, I = 0, phase shifts 5 are given as 
a function of the energy in form of the effective range functions. For n — d 
scattering this function is defined as K{E ) = k cot 5, with E° = E — E^. The 
solid line in the figures represents this function in the interval [0, l-E^I]- The 
solid points in the figures are the results obtained from the integral relations 
using bound state like wave functions at the corresponding energies. As can 
be observed, the results using the bound state like wave functions are in com- 
plete agreement with the exact results given by the solid line in all the energy 
interval. 
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Figure 1. The n — d effective range function for J = l/2 + (left panel) and J = 3/2 + (right 
panel) . The solid points are obtained from the second order estimates using the integral relations 
with bound state like wave functions at the corresponding energies. 

3 Analysis of Three Nucleon Force Models 

In order to reproduce correctly the three-nucleon bound state energy, differ- 
ent three-nucleon force (TNF) models have been constructed during the past 
years as the Tucson-Melbourne (TM), Brazil (BR) and the UrbanalX (URIX) 
models \12\ [TBI [T4"] . These models are based on the exchange mechanism of 
two pions between three nucleons. More recently, TNFs have been derived [15J 
using a chiral effective field theory at next-to-next-to-leading order. A local 
version of these interactions (hereafter referred as N2LOL) can be found in 
Ref. [16J. At next-to-next-to- leading order, the TNF has two unknown con- 
stants that have to be determined. It is a common practice to determine 
these parameters from the three- and four- nucleon binding energies (-B( 3 H) 
and B( 4 He), respectively). 

The n — d doublet scattering length 2 a n d is correlated, to some extent, 
to the A = 3 binding energy through the so-called Phillips line [TTJ, ITS] . 
However the presence of TNFs could break this correlation. Therefore 2 a n d 
can be used as an independent observable to evaluate the capability of the 
interaction models to describe the low energy region. In Ref. [19] results for 
different combinations of NN interactions plus TNF models are given. These 
results are shown for the quantities of interest in Table I and are compared to 
the experimental values of the binding energies and 2 a n d [20J. From the table, 
we can observe that the models are not able to describe simultaneously the 
A = 3, 4 binding energies and 2 a n d- 

In Ref. [3] a comparative study of the aforementioned TNF models has 
been performed. Let us briefly review their structure. From the general form 



W= W(i,j,k) 

i<j<k 



(7) 
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Table 1. The triton and 4 He binding energies B (MeV), and doublet scattering length 2 <w (fm) 
calculated using the AV18 and the N3LO-Idaho two-nucleon potentials, and the AV18+URIX, 
AV18+TM' and N3LO-Idaho+N2LOL two- and three-nucleon interactions. The experimental 
values are given in the last row. 



Potential 


5( 3 H) 


5( 4 He) 




AV18 


7.624 


24.22 


1.258 


N3LO-Idaho 


7.854 


25.38 


1.100 


AV18+TM' 


8.440 


28.31 


0.623 


AV18+URIX 


8.479 


28.48 


0.578 


N3LO-Idaho+N2LOL 


8.474 


28.37 


0.675 


Exp. 


8.482 


28.30 


0.645±0.003±0.007 



a generic term can be decomposed as 

W(l, 2, 3) = aW a (l, 2, 3)+bW b (l, 2, 3)+dW d (l, 2, 3)+c D W D (l, 2, 3)+c E W E (l, 2, 

(8) 

Each term corresponds to a different mechanism and has a different operato- 
rial structure. The specific form of these terms in configuration space is: 

W a (l, 2, 3) = Wq(ti ■ T 2 )(eri • r 3 i)(cr 2 • r 23 )y(r 3 i)y(r 23 ) 
W 6 (l,2,3) = Wo(ri • T 2 )[((T 1 ■ a 2 )y(r 31 )y(r 23 ) 

+ (<Ti • r 3 i)(<r 2 • r 23 )(r 3 i • r 23 )t(r 3 i)t(r 23 ) 

+ {tr\ ■ r 3 i)(<r 2 • r 31 )t(r 31 )y(r 23 ) 

+ Ol • r 23 )(a 2 ■ r 23 )y(r 31 )t(r 23 )} (9) 
W d (l,2,3) = W (t 3 ■ n x r 2 )[(cr 3 • a 2 x cri)y{r 3 i)y{r 23 ) 

+ (o-i • r 3 i)(<T 2 • r 23 )(cr 3 • r 3 i x r 23 )t(r 31 )t(r 23 ) 
+ (<Ti • r 31 )(a 2 ■ r 31 x <r 3 )t(r 31 )y(r 23 ) 
+ (cr 2 • r 23 )(cr 3 • r 23 x o-i)y(r 3 i)t(r 23 )] , 

with Wo an overall strength. The b- and d-terms are present in the three 
models whereas the a-term is present in the TM' and N2LOL and not in 
URIX. The last two terms in Eq.([8]) correspond to a two-nucleon (2N) contact 
term with a pion emitted or absorbed (D-term) and to a three-nucleon (3N) 
contact interaction (E-texm). Their local form, derived in Ref. |16j . is 

W D (1,2,3) = ^<f (n • T 2 ){( < r 1 • a 2 )[y(r 31 )Z (r 23 ) + Z (r 31 )y(r 23 )] 
+ (cri • r 31 )(a 2 • r 3 i)t(r 3 i)Z (r 23 ) 
+ (o"l • ^23)(o-2 • r 2 3)^o(^3l)i(^23)} 

We(1,2,3) = W e {t 1 ■ T 2 )Z (r 31 )Z (r 23 ) . 

The constants W D and Wq fix the strength of these terms. In the case of the 
URIX model the D-term is absent whereas the i?-term is present without the 
isospin operatorial structure and it has been included as purely phenomeno- 
logical, without justifying its form from a particular exchange mechanism. The 
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Figure 2. The p-d analyzing powers A y and iT n at S p = 3 MeV for the AV18+N2LOL (grey 
band) and AV18+UR (solid line) models. Experimental data are from Ref. |21| . 

different form for the profile functions of each model, y(r), t(r) and Zq(t) are 
given in Ref. [3] . In that reference the strengths relative to the different terms 
have been varied in order to reproduce, as close as possible, 5( 3 H) and I?( 4 He) 
and 2 a n d- With these new parametrizations selected polarization observables 
can been calculated and compared to the experimental data. In particular 
using the N2LOL three-nucleon force a small improvement in the description 
of the vector analyzing powers A y and iT\\ at low energies has been obtained. 
This is shown in Fig. [2] in which the predictions for p — d A y and d — p iTxx at 
E p = 3 MeV of the AV18+N2LOL model (grey band) is compared to those 
ones of the AV18+UR model (solid line). From the figure we can observed 
that the discrepancy has been appreciable reduced. 

The effects of the N2LOL TNF is much more evident in the four nucleon 
system. In fact, in Fig. [3]the p— 3 He analyzing power A y is shown at three en- 
ergies using the N3LO-Idaho NN interaction plus the N2LOL TNF (dashed 
line) and the AV18+UR model (solid line). We can see the big effect pro- 
duced by the inclusion of the N2LOL TNF model. It should be noticed that 
the observable calculated using the N3LO-Idaho or the AV18 NN forces alone 
results close to the predictions of the AV18+UR model (see Ref. [22]), indi- 
cating that the improvement is given by the inclusion of the N2LOL force. 
As in the three-nucleon system, this TNF model considerable reduces the 
discrepancy obtained in the description of this observable. 

4 Conclusions 

Two different aspects of the description of few-nucleon systems have been 
discussed. Firstly, scattering states below the deuteron breakup threshold has 
been calculated using bound state like wave functions. The starting point 
in this analysis was the integral relations recently derived from the KVP. 
Finally an analysis of the effects of TNF models has been briefly discussed 





Figure 3. The p— 3 He analyzing powers A y at three energies for the N3LO-Idaho+N2LOL 
(dashed line) and AV18+UR (solid line) models. Experimental data are from Refs. [231124] . 



in the vector analyzing powers in p — d and p— He scattering. In particular 
it was shown that the inclusion of the N2LOL TNF appreciable improves 
the description of those observables. Further studies on these subjects are at 
present in progress. 
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